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Abstract 


A fast computer program has been developea that 
can be used in two basic modes: (1) an analysis 

mode for steady, transonic, potential flow through a 
given planar cascade of airfoils and (2) a design 
mode for converting a given cascade into a shockless 
transonic cascade. The design mode can automati- 
cally be followed by the analysis mode, thus con- 
firming that the new flow field fuund is shock 
free. The program generates its own multilevel 
boundary-conf orming computational grids and solves a 
full-potential equation in a fully conservative 
form. The shockless design is performed by imple- 
menting Sobieciky's fictitious-gas elliptic con- 
tinuation concept. 

Nomenclature 

a speed of sound (isentropic) 

a* speed of sound (critical 

if speed of sound (fictitious) 

Cq coefficient of aerodynamic drag force 

(x direction) 

C|_ coefficient of aerodynamic lift force 

(y direction) 

c airfoil chord lengtn 

D detenninant: *(*,y)/8(X,Y) 

g y-distance between corresponding points 

on neighboring airfoils 
M local Mach number (M ■ q/a) 

M* critical Mach number (M* » q/a») 

Mj Mach number at upstream infinity 

M^ Mach number at downstream infinity 

m coordinate direction orthogonal to 

streamline 

P constant in fictitious-gas relation 

q magnitude of local velocity vector 

S entropy 

U,V contravariant components of velocity 

vector in (X,Y) plane 

u.v components of velocity vector in (x,y) 

pi me 

x,y Cartesiar coordinates in physical plane 

X,Y Cartesian coordinates in computational 

plane 

ai.ap free-stream angles at upstream ano down- 
stream infinity 

e cascade stagger angle 

f ratio of specific heats 

« angle between x axis ano velocity vector 

V Prandt 1-Meyer function 


p isentropic fluid density 

0 * critical fluid density 

pf fictitious fluid density 

# velocity vector potential (T(> • 7) 

stream functioii 

Introduction 

In the general case of transonic cascade flow, 
supersonic regions terminate with shocks. These 
shocks create vorticity and generate entropy in a 
flow field that was initially irrotational and 
homentropic. As a consequence the aerodynamic drag 
force sharply increases (wave drag) and the total 
energy decreases, resulting in a rapid decay of the 
aerodynamic efficiency of the cascade and an abrupt 
increase in the aerodynamic noise level. In many 
experiments it has been observed that, if the Mach 
number just ahead of the foot of the shock wave is 
larger than approximately 1.3, the boundary layer 
starts to separate, leading to complex and poten- 
tially dangerous unsteady flow phenomena and mech- 
anical vibrations. 

Choked flow represents yet another undesirable 
phenomenon associated with transonic cascade flow. 
Choking places an upper limit on the mass flow 
through a given cascade. As a countermeasure the 
airfoils in the cascade are often positioned farther 
apart, decreasing cascade solidity. This results in 
? decrease in flow turning angle through the cascade 
ond a drop in pressure rise across the cascade. 

The main objective of this work is therefore to 
eliminate the shocks (and possibly even the chokeo 
flew) by slightly altering portions of the contour 
of a given airfoil ir the cascade. 

Analysis 


Governing Equations 

This work is based on the fictitious-gas con- 
cept of Sobieczkyl and the full-potential, steady, 
transonic turbomachinery analysis codes of 
Dul ikrav ich. 2 The analysis was derived exten- 
sively in earlier works^»^ and will be repeated 
here in its concise form only. 


In the case of a steady, two-dimensional, 
irrotational isentropic flow of an inviscid, com- 
pressible fluid the conservative form of the con- 
tinuity equation is 
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(pu) + (ov) ■ 0 

I * •/ 

Equation (1) can also be expressed in its non- 
conservative full-potential form 


( 1 ) 
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( 2 ) 


or in its canonical operator formS.2 

^p(l - - p(l - ■ ° 


(3) 


Equation (3) represents a quasi-1 inear, second- 
order partial differential equation of mixed 
elliptic-hyperbolic type that accepts isentropic 
discontinuities in its solution. These isentropic 
shocks satisfy mass conservation 




(4) 


and differ from the Rank ine-nugoniot shocks 
(Table 1). Superscript H in equation (3) desig- 
nates upstream differencing, and superscript E 
designates central diffe-encing to be used for the 
evaluation of particular second derivatives. Solu- 
tion of this steady-state equation is obtained as an 
asymptotic solution to an artificially unsteady^ 
equat ion 

o((l-M^)0 S, *0,^*^Ca^St*^''^mt*‘^t)■O 


Consequently the fully conservative form of the con- 
tinuity equation (eq. (1)) becomes 

J^(pUD * «x' ^ J ° 

where the artificial viscosity terms «x«*Y 
represent principal parts of a truncation erro'^ of 
equation (3). 

The computational grid in the x.y plane is 
generated by usinq a sequence of simple geometric 
transformations^'lO incorporating a single 
conformal-mapping function, elliptic polar coordi- 
nates, and nonorthogonal coordinate stretching and 
shearing. The uniform grid (Fig. 2) in the computa- 
tional (X,Y) plane thus remaps back into the body- 
fitted, quasi-orthogonal grip of figure 3 in the 
physical (x,y) , »ne. The iterative solution pro- 
cess of equation (3) is accelerated by using a four- 
level, consecutive-grid refinement procedure. 

All the flow parameters are nondimensionalized 
with respect to the critical conditions denoted by 
an asterisk so that the isentropic relations used 
for t:.e local fluid density and the speed of sound 

.rell 
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( 10 ) 


(5) 

for large times, where E, n, and c are coeffi- 
cients. This equation is solved by using an itera- 
tive line overrelaxation whe'^e consecutive iteration 
sweeps through the flow field are considered as 
steps in an artificial time direction The steady 
part of the residual (or error) of equation (b) is 
..Iways evaluated by using equation (1) supplemented 
by a directional numerical viscosity in a continu- 
ously fully conservative form, thus uniquely captur- 
ing possible isentropic shocks. 

For the purpose of a type-dependent,^ 
rotated^ finite difference evaluation of the 
derivatives in equation (3) and a finite area^ 
evaluation of the first derivatives in equation (1), 
the flow field and the governing equations are 
transformed from the physical (x,y) plane (Fig. 1) 
into a rectangular (X,Y) computational domain (Fig. 

2) by using local isoparametric bilinear mapping 
functions. 

If the geometric transformation matrix is 



then the contravariant velocity components in the 
(X,Y) plane are 

(7) 

V,yJ 


Shock-Free Surface Design 

Within the last decade several versions of an 
indirect (hodograph) design approach based on 
Garabedian's method of complex characteristics have 
been published. 12 j^e method proved to be a 
powerful tool for the design of high-performance 
airfoils and cascades, but handling the complicated 
boundary and initial-vxlue problems in a four- 
dimensional computational space for practically 
interesting design case, requires a large amount of 
experience. It is therefore desirable to develop 
efficient direct - or nearly direct-design methods. 

This task can be accomplished by prescribing a 
smooth, shock-free pressure distribution along a 
portion of a given airfoil contour in a cascade and 
then determinir.g a partially new airfoil shape con- 
sistent with the prescribed surface flow condi- 
tions. Because of the highly nonlinear character of 
the transonic flow this design technique generally 
does not provide an entirely shock-free flow field. 
In order to completely eliminate all shocks (and the 
associated wave drag) from the flow field, a number 
of suen designs must be performed, and an optimiza- 
tion technique must be devised to search for a cas- 
cade that maintains an entirely shock-free flow 
field for a specific set of flow parameters. 

Shock-Free Flow Field Design 

Shock-free, or shock less, flow means that the 
fluid decelerates from a supersonic speed to a sub- 
sonic speed not discontinuously (shocked flow), but 
smoothly over a finite distance (isentrop c recom- 
pression). This requires determination of a modi- 
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fled supersonic region that Is longer but flatter 
and thus can possible unchoke a flow that was 
originally choked. The sonic line bounding such a 
supersonic region in an otherwise subsonic flow 
field must not have inflection points if shocks are 
to be avoided. 

To eliminate the possibility of obtaining 
shocks anywhere in the flow field, Sobiectky pro- 
posed^ and successfully appliedl3,14 the concept 
of a fictitious-gas, shock-free design that corre- 
sponds to an elliptic continuation^ from the sub- 
sonic flow into local supersonic flow domains. This 
design technique uses isenlropic relations for fluid 
density (eg. (9)) and the sound speed (eq. (10)) 
only where the flow is locally subsonic. At every 
point where the flow is locally supersonic, modified 
(artificial) relations are used for the fluid den- 
sity and the sound speed so that the governing equa- 
tion remains elliptic throughout the flow field. 
Therefore any conservative computer code capoble of 
solving a subcritical potential flow field can be 
modified to include the fictitious gas concept and 
then serve as a tool for determining the sonic line 
shape. 

It IS important to point outl^ that the flow 
field outside the supersonic bubble calculated from 
the f ictitious-gas relations is already the correct 
subsonic flow field. It is only the supersonic part 
of the flow field that must be recomputed ana from 
this recalculation a new portion of the shock-free 
airfoil surface determined. Lift ana drag coeffi- 
cients are also already design results, and they 
will not be altered by the subsequent recomputation 
of the local supersonic region. The sole purpose of 
originally using an fictitious gas (nioaified density 
and sound speed relation) is thus to determine a 
shape of the sonic line that is compatible with an 
entirely shocn-free flow field. 



when M > 1. It IS desirable to use a formula for 
pf/p* that includes a single (preferably con- 
stant) input parameter P tha^ makes the fictitious 
gas more or less compressible in the supersonic 
regions. Such a function must not have minimum or 
maximum values in the range of expected relative 
local Mach numbers, because at such points the local 
fictitious speed of sound af/a« is infinite. 

This can be observed if the general continuity 
criterion 



is used to obtain the relation for if la*. 

After taking a logarithm and a derivative of both 
sides of equation (1?) one gets a general expression 
for the speed of sound of the fictitious gas 



for the purpose of guaranteeing an entirely 
shock-free flow field the values for pf/p« must 
always be higher than the values required by the 
parabolicity condition; that is, (pf/p«) • 

(Fig. 5). The final condition tor the relation 
pf/p* « F(M«; P) is that it should be a very 
simple function that will also produce a simple ex- 
pression for af/a». 


In tne two-dimensional case of a cascade of 
airfoils the values of potential f and stream 
function i on the sonic line obtained froni the 
f ici it ious-gas calculation serve as the initial data 
for an integration by the method of characteris- 
tics. This integration is performed in a triangle 
EKC of a rheograph (Fig. 4) plane (local flow angle 
0 versus Prandt 1-Meyer function v(M«). This 
recalculat ion (now using isentropic gas relations 
(eq. (9) and Eq. (10)) is performed only inside the 
supersonic flow domain. The new section of the air- 
foil surface is then determined from the arc i . 0, 
which is interpolated froni the solution in the char- 
actens'ic triangle EMC (Fig. 4). The new airfoil 
is slightly flatter than the original shape: Isen- 

tropic flow requires more space to pass than the 
fictitious one. 


In the present work we use the relation 


Of 1 - V^l ♦ 4P{M, - 1) 

^.1 (14) 


which gives 
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? / 1 - \ 
? ■ ^ V j 


(15) 


The fictitious-gas technique is not limited in 
application only to two-dimensional planar problems 
like a hodograph technique; it can be successfully 
applied to both arbitrary two-dimensional and three- 
dimensionall'^*^^ conf igurat ions. 

Fictitious-Gas Relations 

The fictitious-gas relation pf/p« is ap- 
plied only in the regions where M > 1, An arbi- 
trary analytic expression for pf/p« is never- 
theless subject to several constraints (Fig. 5). It 
should satisfy the first-order continuity condition 
on the sonic line in the flow field; that is. 


Results 

Based on the preceding analysis, computer 
code^^ DCAS2D has been developed and tested for 
the following sequence of test cases. For the pur- 
pose of illustrating basic features of the flow 
through planar cascades of airfoils the flow around 
an isolated NACA 001? airfoil in free air and the 
flow through a cascade of NACA 0012 airfoils were 
analyzed. Airfoils in the cascade had zero stagger 
angle (e • 0*) and a gap-chord ratio of 3.o 
(g/c»3.6). The free-stream angle was zero at both 
upstream and downstream infinity. In the case of an 
incompressible free stream (Mj ■ 0.001) the result 
obtained for the cascade did not differ from the 
result obtained for an isolated airfoil.^ But in 
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the case of a transonic flom (Mj ■ 0.8) the cas- 
cade effects (Fig. 6) were very significant even for 
such widely spaced airfoils. 

To demonstrate the applicability of shock-free, 
fictitious-gas design to realistic lifting, stag- 
gered cascades, we use a simple analytical shape 
generator for geometry definition of the input air- 
foils. Flexible geometry definition is most useful 
for parametric studies of cascades. Here we use a 
formula lor blade section oefinition: 

y . Ax ♦ ♦ x^'^d - x)^'^ (d ♦ Ex ♦ Fx(l - x)) 

(16) 

with a proper choice of the parameters to control 
leading- ood trai 1 mg-edge radii, angles, and thick- 
ness distribution. 

Figure 7 shows a cascade of this family, with a 
sonic line and the corresponding modification on ihe 
upper surface, where shock-free redesignl6 for 
chosen operating conditions flattened the a’-';' 1 
contour. For this cascade the global flow geo- 
metric garometers^were « C.8,^g/c - 0. 
aj B 41 , • 17 , and B • 27.3 . 

It is shown in Figure 8 that the a' ‘flow 
through this cascade with geometry given by equation 
(16) and a gap-chord ratio of 0.8b is not shock-tree 
but contains a very strong shock. For the same 
global flow conditions the design mode of the code 
0CAS20 IS then used to obtain a new shape (Fig. 7). 
This shape is different from the initial airfoil 
only between 3 and 38 percent of chord on the upper 
surface, resulting in a shock-free (Fig. 9) pressure 
distribution. Then the analysis mode of DCA2S0 is 
used to verify that the design flow agreement is 
excellent (Fig. 10). 

This IS the first of a series of examples from 
a parametric cascade airfoil shape study. 1“ 

Although the new airfoil losses shock-free proper- 
ties at off-design conditions (Fig. 11), the result- 
ing shock IS still considerably weaker than a shock 
on the original airfoil. An optimum cascade for a 
range of operating conditions can be obtained by 
combining the f ict it lous-gas design concept with an 
optimization technique. 

As already mentioned, computer code DCAS2D is 
capable of converting a choked, shocked cascade flow 
field into an uncooked, shock-free flow field. 

To illustrate this feature we selected a non- 
staggered cascade of NACA 0018 airfoils having a 
gai»-chord ratio g/c of 1. Note that a simple one- 
dimensioialU flow assumption predicts that 
the flow through thi-. cascade will choke if 
Hj > 0.577. Therefore we used the design mode 
of the DfASPD code with Mj . 0.582 and the 
fictitious-gas parameter P • 500. The resulting 
flow field (Fig. 12) is unchokeJ and entirely shock 
free. 

All the calculations were performed without 
taking into account viscous boundary layer effects. 
For this purpose one may use a standard boundary 
layer calculation procedure because shock - boundary 
layer interaction effects do not exist in a shock- 
free flow. The V iscous/ inv ISC Id calculation can be 
performed iteratively with a treatment of trailing- 


edge viscous interaction, as has been demon- 
strated” for isolated supercritical airfoils. 

Concluding Remarks 

An efficient and reliable computer program, 
DCAS20, has been developed and tested that automati- 
cally performs partial redesign of a given airfoil 
shape in the cascade for the purpose of eliminating 
shock waves and the associated wave drag. The code 
represents an application of already known and suc- 
cessfully applied numerical techniques for transonic 
flow analysis and the shock-free flow field design. 
These techniques are based on the finite volume and 
a fictitious-gas approach, respectively. A new for- 
mula for the fictitious-gas relation, accompanied 
with the related physical constraints, has been 
suggested. 

The computer code is entirely self-sufficient 
in generating its own multilevel boundary-conforming 
grids. The code can operate separately as a shock- 
free cascade design code and also as a general tran- 
sonic cascade analysis program with the capability 
tn accurately capture isentropic shocks. 
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TABLE I. RANKINE-HUGONICT (Mb*) AND ISENTROPIC 


(Mb*. Mb) SHOCK JUMP RELATIONS 


Ma* 

Mb* 

Mb* 

Ma 

Mb 

1.00500 

0.99502 

8 . <9500 

1.00601 

0 . 9940 , 

l.OlOi 0 

0.99010 

0.99001 

1.01204 

0.98804 

1.01500 

0.98522 

0.98502 

1.01808 

0.98210 

1.02000 

0.98059 

0.98005 

1.02415 

0.97617 

1.02500 

0.97561 

0.97504 

1.05025 

0.97027 

1.05000 

0.97087 

0.97006 

1.05655 

0.96459 

1.05500 

0.96618 

0.96508 

1.04245 

0.95855 

1.04000 

0.96154 

0.96011 

1.04859 

0.95269 

1.04500 

0.95694 

0.95514 

1.05475 

0.94687 

1.05000 

0.95258 

0.95017 

1.06095 

0.94107 

1.05500 

0.94787 

0.94520 

1.06715 

0.95529 

1.06000 

0.94 540 

0.94024 

1.07555 

0.92952 

1.06500 

0.95897 

0.95528 

1.07959 

0.92578 

1.07000 

0.95458 

0.95055 

1.08585 

0.91806 

1.07500 

0 . 95025 

0.92558 

1.09215 

0.91256 

1.08000 

0.92595 

0.92045 

1.09845 

0.90668 

1.08500 

0.92166 

0.91548 

1.10476 

0.90101 

1.04000 

0.91745 

0.91054 

1. 11110 

0.89557 

1 . 045 C 0 

0 91524 

0 90560 

1.11747 

0.88974 

1. 10000 

0.90909 

0.90067 

1.12585 

0.88415 

1.10500 

0.90498 

0.89574 

1.15026 

0.87854 

1.11000 

0.90090 

0.89081 

1.15670 

0.87297 

1.11500 

0.89686 

0.88589 

1.14515 

0.86741 

1.12000 

0.89286 

0.88096 

1.14965 

0.86188 

1.12500 

0.88889 

0.87605 

1.15615 

0.85656 

1 . 1 5000 

0.88496 

0.87115 

1.16265 

0.85086 

1.15500 

0.88106 

0.86622 

1.16920 

0.84558 

1.14000 

0.87719 

0.86152 

1.17577 

0.85991 

1.14500 

0.87556 

0.85641 

1.18256 

0 . 8541.6 

1.15000 

0.86957 

0.85151 

1 . 18898 

0.82905 

1.15500 

0 86580 

0.84662 

1 . 1 V 565 

0.82561 

1.16000 

0.86207 

0.84175 

1.20229 

0.81821 

1 . 16500 

0 85857 

0.85684 

1.20899 

0.81285 

1 . 17000 

0.85470 

0 85195 

1.21571 

0.80747 

1.17500 

0 35106 

0.82707 

1.22245 

0.80212 

1.18000 

0.84746 

0.82219 

1.22922 

0.79678 

1 . 18500 

0 84588 

0 81752 

1.25602 

0.79147 

1.19000 

0.84054 

0.81245 

1.24284 

0 . 78617 

1.19500 

0.85682 

0.80759 

1.24969 

0.78088 

1 . 20000 

0.85555 

0.80272 

1.25656 

0.77561 

1.20500 

0 82988 

0.79787 

1 .26546 

0.77056 

1.21000 

0.82645 

0.79501 

1.27059 

0.76512 

1 .21500 

0.82505 

0 . 78816 

1.27755 

0.75990 

1.22000 

C 81967 

0.78551 

1 28454 

0.75469 

1 . 22500 

0 81655 

0.77847 

1.29155 

0.74950 

1.23000 

0 81501 

0 77565 

1.29840 

0.74455 

1.25500 

0 80972 

0.76880 

1 . 50547 

0.75917 

1 .24000 

0.80645 

0 . 76597 

1.51257 

0.75402 

1 ,24500 

0 80521 

0 .75914 

1.51970 

0.72889 

1.25000 

0.80000 

0.75452 

1.52686 

0 72577 

1.25500 

0.79681 

0.74950 

1 . 55405 

0.71867 

1.26000 

0.7 9565 
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Figu't L - Piinir cii:i0t of ainoiii In phys- 
ical U.yi plant. 
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figurt 2. ■ Ptrioflic flo* fitlfl in (X, VI compuUlionii plant. 



Figurt 3. - Computalional grid in physical (t, yl plant. 
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Fiqur* 6. - Trjnsonic c«sc»0« iflecU. lAirfoil, 
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